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We present a new calculation of the D — > it and D — > K form factors from 



QCD light-cone sum rules. The MS scheme for the c-quark mass is used 
and the input parameters are updated. The results are /^_(0) = 0.67^0 07, 
f+ K (0) = 0.75±g;j£ and /o w (0)//5 Ji -(0) = 0.88 ± 0.05. Combining the cal- 
culated form factors with the latest CLEO data, we obtain \ V c d\ = 0.225 ± 
0.005 ±0.003 tomt and | V^ rf |/| V^ s | = 0.236 ± 0.006 ± 0.003 ± 0.013 where the 
first and second errors are of experimental origin and the third error is due 
to the estimated uncertainties of our calculation. We also evaluate the form 
factors fp n and f^ K and predict the slope parameters at q 2 = 0. Further- 
more, calculating the form factors from the sum rules at q 2 < 0, we fit them 
to various parameterizations. After analytic continuation, the shape of the 
D —* 7T, K form factors in the whole semileptonic region is reproduced, in a 
good agreement with experiment. 



1. Introduction 



Recent measurements of the semileptonic D — » irivi and D —* lf£z^ decays by CLEO 
collaboration |1] [2| provide new accurate results on branching fractions and differential 
decay rates, in addition to the previously accumulated data [HJ HI El [6] . The decay rate 
distributions in the bins of the variable q 2 (the invariant mass squared of the lepton pair) , 
yield the products of transition form factors and CKM matrix elements, | V c d\f^ (q 2 ) and 
\V cs \f^ K (q 2 ). In addition, the form factor shapes are reconstructed and fitted to various 
parameterizations. With the new data available, it is timely to update the theoretical 
analysis of the D — > tt and D — > K form factors, aiming at more accurate determination 
of \V cd \ and \V CS \. 

In this paper, we recalculate the D — > tt,K form factors from QCD light-cone sum rules 
(LCSR's). In this method [TJ, the correlation function of quark currents is constructed 
in a form of a transition matrix element between the vacuum and the final hadron 
state. In our case, the quark current with D-meson quantum numbers is correlated 
with the charmed weak current, whereas the hadron state is the on-shell pion or kaon. 
Two different representations of the correlation function are then equated: the operator- 
product expansion (OPE) near the light-cone and the dispersion integral over hadronic 
states. In the latter, the ground D-meson contribution containing the D — > tt or D — > 
K form factor is singled out. Applying the quark-hadron duality approximation, the 
remaining dispersion integral over the higher states is replaced by the integral over 
the OPE spectral density. The LCSR approach, though having a limited accuracy, 
provides analytical expressions for the form factors, in terms of finite quark masses and 
universal light-cone distribution amplitudes (DA's) of pion or kaon. Importantly, the 
heavy-light form factors calculated from LCSR's with gluon radiative corrections, include 
both "hard-scattering" and "soft-overlap" components, and the latter is predicted to be 
the dominant one. 

The fact that the correlation function is calculated at a finite heavy quark mass, simpli- 
fies our task, because the OPE of 6-quark and c-quark currents have the same analytical 
expressions. Only the quark mass value has to be changed and the normalization scales 
have to be adjusted. Hence, for example, the LCSR for fp n represents a by-product of 
the LCSR obtained for . The latter form factor used to determine \ V u b\ is the most 
familiar application El Q33 EH H3 H3] of this method 0. In what follows, we employ 
the recent update of LCSR for B —* ir form factors presented in [13J . Importantly, the 
D — > tt, K form factors obtained from the sum rules, being confronted with experimental 
data, provide an important test of the whole method. 

Compared with the previous calculations of D — > tt,K form factors from LCSR's 
|15l I16j . certain modifications and improvements are made. First of all, following |13j . 
we systematically use the MS scheme for the c-quark mass, whereas earlier calculations 
switched to the pole mass in the final sum rules. In this respect we benefit from recent 
accurate determinations [T7] of the c-quark mass from the charmonium QCD sum rules. 



1 Interestingly, one of the earliest applications [H] of the "sum rules on the light-cone" was to the 
D —* form factors. 
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Also the strange quark mass entering the sum rule for D — > K form factor has a smaller 
uncertainty than before. Furthermore, we use the improved determination [18] of the 
SU (3) fi violating Gegenbauer coefficient of the twist- 2 kaon DA, and the updates of the 
pion and kaon twist-3,4 DA's from |19j . 

The main novelty in this paper concerns the q 2 -dependence of the form factors. First, 
we determine the slope parameters at q 2 = 0, which involve the second form factor 
^Dtt(K)' l a tter is calculated using the same method and input. Note that for 

D — > tt(K) transitions, LCSR's are applicable only in the lower part of the region 
< q 2 < (rriD — m n ^ K ^) 2 , accessible in the semileptonic D — > ir(K)lvi decays. At 
q 2 approaching m 2 , the virtuality of the c-quark becomes a soft scale, and OPE is 
not reliable. In this paper we predict the form factor shapes, combining LCSR with 
analyticity. We employ the spacelike region q 2 < 0, where the light-cone OPE works 
even better, than at small positive q 2 . The LCSR results for the form factors at q 2 < 
are fitted to various parameterizations, such as the dispersion relation with an effective 
pole |20j and the recent version of series parameterization [21] . We then make use of the 
analytic continuation from negative to positive q 2 and predict the form factors in the 
whole semileptonic region. 

The paper is organized as follows. In Section 2 we present the outline of the method 
and discuss the expected accuracy of our calculation. In Section 3 the input parameters 
are listed and in Section 4 the numerical results for the form factors at zero momentum 
transfer are presented. In Section 5 we compare the calculated form factors with exper- 
imental data, and determine \V CC [\ and \V CS \. In Section 6 we turn to the determination 
of the form factor shape and present the fit results for various parameterizations. The 
calculated shapes are used to predict the total semileptonic widths. Finally, Section 7 
contains a concluding discussion. In Appendix A the relevant definitions of light-cone 
DA's and their expressions used in LCSR's are collected, and in Appendix B the formulae 
for the contributions to LCSR's are presented. 

2. Outline of the LCSR method 

The central objects of our interest are the D — > tt form factors f^ w (q 2 ) defined in a 
standard way from the hadronic matrix element: 

<7T-(p) | d^c | D°(p + q)) = 2/+ r (< ? 2 )p / , + (f^(q 2 ) + /rf))^ , (1) 

and the analogous form factors fp K (q 2 ) of D° — ► K~ transition, obtained by replacing 
d — > s and tt~ — * K~ in the above. In what follows, we work in the isospin symmetry 
limit, so that the D + — » K° and D° — > K~ hadronic matrix elements are equal and the 
D + — ► 7T° form factors are obtained by multiplying fp n (q 2 ) with l/y/2. The form factor 
/~ is combined with / + in the scalar form factor 

f D « { K)tf) = fL ( K)(<l 2 ) + m 2 ^2 W)^ ' (2) 

m D m -K(K) 

which plays a minor role in semileptonic transitions and is "visible" only in D — » Kfiu^. 
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The correlation function used to derive LCSR's for D —* it form factors is defined as 

Fffaq) = iy"d 4 xe^^(p)|T{ ( i>)7^ C ( a ;),m c c(0)i7 5U (0)}|0) 

= F 7r (q 2 ,(p + q) 2 )p fl + F^(q 2 ,(p + q) 2 )q f , . (3) 

Replacing d — ► s and tt — > K in the above, we obtain , the correlation function for 
D ^ K form factors. The two invariant amplitudes F'*^ and i* 1 *^) yield two separate 
sum rules for and (f^K) + /Djr(iC))' respectively. Using Q, one then obtains 

/dttW so that in our calculation, by default, /^(^ (0) = Z^^C )- 

At <C m 2 and (p+g) 2 ^ m e> the c-quark propagating in the correlation function has 
a large virtuality and the product of the c-quark fields is expanded near the light-cone 
x 2 ~ 0. This expansion starts in leading order (LO) from the free c-quark propagator 
and includes the 0(a s ) (NLO) corrections due to hard-gluon exchanges between the 
quark lines and soft-gluon emission. A more detailed derivation can be found e.g. in [9j, 
where also the origin of the light-cone expansion in the correlation function is explained. 

The OPE result for the correlation function ([3]) is cast in a factorized form, where 
the perturbatively calculable kernels are convoluted with the pion DA's of growing twist 
t = 2, 3, 4, ... For the invariant amplitude F n one obtains: 

(p + <!?)] QpE = £ F ^(q 2 , (p + qf) + ^ £ Ff(q 2 , (p + q) 2 ) , (4) 
t=2,3,4 i=2,3 

where the LO (NLO) parts in a s are the convolutions 

F^tfAp + qf) = J %r ( J ) ( g 2 ,(j)+g) 2 ) m^ ) M)^(«,M)- (5) 

The perturbative kernels T^l stem from the c-quark propagator, and include the 

loops with hard-gluon exchanges in 0(a s ). The pion DA's <j>$ (u, fj) represent universal 
vacuum-pion matrix elements of light-quark and gluon operators. The simplest bifocal 
quark-antiquark operators d(x)T a u(0) (where T a is a generic combination of 7-matrices) 
originate after contracting the free c-quark fields in the correlation function (J3|) . In 
addition, soft gluons emitted from the propagating c-quark, together with light quarks 
and antiquarks, form DA's of higher multiplicity, starting from the three-particle (quark- 
antiquark-gluon) DA's of t = 3,4. In ([5]), the integration over u is a generic notation for 
the momentum distribution between constituents of the pion in DA's. The definitions 
and explicit expressions for all relevant DA's are given in Appendix A. The presence 
of the terms with 2 < t < 4 in @) reflects the currently achieved accuracy in OPE: 
the twist 2,3,4 terms in LO (including three-particle contributions of twist 3,4) and the 
twist- 2 and twist-3 (two-particle DA's) terms in NLO. 

The factorization scale [i in ([5]) separates the perturbative kernels dominated by near 
light-cone distances from the long-distance quark-gluon dynamics in DA's. The collinear 
divergences in the twist- 2 of OPE are absorbed in the logarithmic evolution of DA's as 
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shown in \10\ lllj . The factorization for the twist-3 part is proved for the asymptotic 
DA's in [12] and confirmed in [13] . 

Note that the two terms in the c-quark propagator proportional to m c and j> c yield two 
different parts of OPE with even (t = 2,4, ..) and odd (t = 3,5, ..) twists, respectively, 
corresponding to different chiralities of r a matrices in the light-quark operators. Hence, 
the twist expansion goes in even and odd twists separately. The two most important 
LO contributions of twist- 2 and twist-3 two-particle DA's to the correlation function 
originate from the ~ m c and ~ p" c parts of the free c-quark propagator, respectively, and 
have simple expressions: 



j-iTT.SHwopart. 
^0 



m 2 — q 2 u — (p + q) 2 u 
o 

i 



\q 2 ,{p + qf) = Un^rric I — j^— — -^{^(u) 

/ mj. — q z u — (p + q) z u I 

2 2 

(6) 



2 + m c + 9 



6 L m 2 — q 2 u — (p + q) 2 u 

where fj, n = m 2 /(m u + ma) and u = 1 — u. In the above, ip^ and c/)^ n , 4>^ v are the 
pion twist- 2 and twist-3 two-particle DA's, respectively. Note that the formal l/m c 
suppression of Fq' 3 versus Fq' 2 is overwhelmed numerically, because the enhanced light- 
quark parameter fj, n is larger than m c . The remaining twist-3 term in LO, due to the 
three-particle DA, not shown in ([6]), is strongly suppressed by the ratio of the small 
normalization factor fy n to m c . 

Turning to higher- twist (t > 3) contributions, one has to take into account that in 
the light-cone OPE each two units of twist are accompanied by an extra x 2 , yielding an 
additional power of the denominator 

D = —o T^~< n5-i (° <«<!)> (7) 

m z — q z u — [p + q) z u 

in the correlation function. Hence, for example, the contribution Fq' 4 of the twist-4 
two- and three-particle DA's is subleading with respect to the twist-2 part F^' 2 in Q, 
with a suppression factor ~ 8 2 (D) , where (D) is the weighted (over DA's) average of 
the denominator ([7]), and 5 2 ~ A.q CD is the normalization factor of the twist-4 DA's. 
Similarly, the twist-5 contributions, not yet included in the currently used version (jlj) of 
OPE, are expected to be suppressed with respect to the twist-3 terms, parametrically 
~ Aq CD (D) . We emphasize that the suppression of higher- twist contributions is effective 
only if both external momenta squared q 2 and (p + q) 2 in ([7]) are kept <C m 2 , that is, 
when the c-quark is sufficiently virtual. 

The detailed calculation of the correlation function, including 0(a s ) corrections to 
the twist-2 and twist-3 part, is given in previous papers and we will not repeat it. Here 
we use the recent update in terms of MS heavy quark mass presented in [131 . where nib 
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has to be replaced by m c , with a corresponding adjustment of the scales. As in [13], we 
use a universal normalization scale /i for a s , quark masses and DA's . The u-, d-quark 
masses and m 2 are neglected everywhere, except in the parameter fi^. 

The correlation function for the D — > K form factors includes SU(3)fi violation 
effects of different origin, starting from 0(m s ) ~ O(m^). One effect is purely kinemati- 
cal, due to the presence of p 2 = m 2 K ^ in the c-quark propagator. The expressions for 
the perturbative kernels have to be modified, so that the term m 2 K uu has to be added 
to the denominator ([7]). Furthermore, there are SU(3)fi violating corrections in DA's. 
In twist 2, the deviations of (fx from <p n is due to differences between the normalization 
parameters (fx vs f n ) and Gegenbauer moments. Importantly, one has to take into 
account the first Gegenbauer moment of the kaon DA, ^ 0, responsible for the mo- 
mentum distribution asymmetry between the strange quark and nonstrange antiquark 
in the kaon. In general, also af ^ af- ^ n twist 3 and 4 DA's one has to replace the 
normalization parameters: fi w — > fife = m K /(m u + m s ), — > fsx and 5 2 — > 8 K , re- 
spectively. In addition, there are 0{m 2 K ) admixtures of the twist-2 DA in the kaon DA's 
of twist 3 and 4. All these and some other less important SU (3) fi violating effects are 
included in the expressions for the kaon DA's [H] presented in Appendix A. 

To access the D — > ir, K form factors, one equates the OPE result to the hadronic 
dispersion relation for the correlation function in the variable (jp + q) 2 , the momentum 
squared of the D-meson interpolating current. In the D — > tt case, the resulting relation 
for the invariant amplitude F 77 is: 



oc 



lmF n (q 2 ,s) 



F (q ,(p + q) ) QpE - m 2 D _ {p + q? +-J ds s _ {p + q? ■ (8) 



h 



In the above, the Z}-meson pole term (with the mass mo and the decay constant defined 
as (0| m c g«75C \D) = m 2 D fr)) : contains the desired form factor f£ n . Also in (jHJ) the 
quark-hadron duality approximation is applied, replacing the hadronic spectral density 
of the higher states by the OPE spectral density. The latter approximation introduces 
the threshold parameter Sq which is non-universal and has to be determined for each 
sum rule independently. Note that Sq is an effective parameter, not necessarily equal to 
the lowest hadronic continuum threshold (mo* + ^7r) 2 - 

After subtracting the integral on r.h.s. of (|Hj) from both sides of this equation, one 
performs Borel transformation of dH) replacing the variable (p + q) 2 with the Borel pa- 
rameter M 2 and exponentiating the denominators. E.g., the powers of the denominator 
D in ([7]) transform as: 

1 2 2- 

W-W^} = ( w - l)! M n( M 2)n-i e • (9) 

Finally, the LCSR for the D — * n form factor is obtained: 

m 2 /A/ 2 / n \ 

Zm DJD \ t= n 3 A ^ t=2f3 / 
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where the functions i* 1 ^'* (q 2 , M 2 , Sq) and F^'* (q 2 , M 2 , Sq) are derived applying the Borel- 
and-subtraction procedure to each twist component of the OPE in (J2J): 

oo 

F^ } (q 2 , M 2 , s°) = B^^M^F^q 2 , (p + q) 2 )} -\Jds e'^lmF^q 2 , s). 

(11) 

The expressions for Fq' and F^' in ()10p are obtained from the corresponding expressions 
given in [13] for I? — > ir LCSR, replacing b — > c and B ^ D. The second LCSR for 
the combination (/^ + /xL) is obtained from the invariant amplitude F^ and has the 
same form as (|10|) , with the invariant amplitudes > replacing F^^\ , an d without the 
factor 1/2 in the coefficient. 

For D — » X form factors, as explained above, the SU(3)fi violation effects are taken 
into account in both LCSR's for f^ K and (/^^ + i n the LO part, keeping p 2 = 

m 2 ^ ^ in the hard kernels and taking into account the 0{m s ) ~ 0{m 2 K ) effects 
in the kaon DA's. For the sake of completeness, we do not expand these expressions 
in m 2 K ~ m s , although it is clear that only the first-order terms of this expansion 
are important numerically. Furthermore, having in mind that both 0(a s ) and 0(m s ) 
corrections are reasonably small, we do not take into account the combined 0{a s m s ) 
effects, originating from nonzero m s and p 2 = m 2 K in the NLO diagrams. These effects 
demand a dedicated calculation, taking into account the mixing between various twist 
components at the 0(m s ) level. Hence, in the adopted approximation F^' 2 ^ = f^' 2 ^ 3 \ 

In Appendix B, we present the expressions for all LO terms j?^'* -2 ' 3 ' 4 entering LCSR's 
for D — > K form factors. In the limit m s — ► {m 2 K — > 0) and fix — * ^ the corresponding 
terms i?^'* -2,3 ' 4 i n the D — > ir LCSR's are easily reproduced. The expressions for the 
NLO terms of twist- 2 and twist-3 are presented in [T3] in a form of dispersion integrals, 
and we replace m& — ► m c . These very bulky formulae are not displayed here. 

The LCSR's for D — > K form factors were also compared to the recent update of the 
LCSR's for B — > K form factors presented in [22j . The expressions presented in Ap- 
pendix B agree with those in [22], up to the twist-3 three-particle part of LCSR, which is 
incomplete in [22J. We also modified the m^-corrections to the two-particle twist-4 DA's 
(see discussion in Appendix A), hence, there are small differences in the corresponding 
twist-4 terms in LCSR's. Both differences have a minor impact on numerical results. 
Finally, we compared our expressions with the LCSR for the D — > P (P = ir, K) form 
factor in LO at q 2 = presented in [16] , where only the 0{m 2 p ) terms are retained. In 
the limit m 2 P = 0, the equations ([731) . (1761) . and ([77)) are reduced to the corresponding 
terms in Eq. (2) of [16| . except we obtain an opposite sign of the contribution of twist- 
4 DA ip4 ; p. We also were not able to completely reproduce the 0{m 2 P ) terms in this 
equation. 

The expected accuracy of LCSR fjlOf) is determined by the uncertainties of the com- 
bined expansion in twists and a s of the correlation function. The choice of Borel param- 
eter plays an important role here. Note that Borel transformation effectively replaces 
the powers of the denominator D in the higher-twist contributions by the inverse pow- 
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ers of M 2 , as seen from Q p. Hence, although M 2 is an arbitrary scale, it should be 
taken sufficiently large to preserve the twist hierarchy. Parametrically, in the limit of 
the heavy c-quark mass, the relevant scaling relation is M 2 ~ 2m c r. Importantly, also 
the scale r has to be much larger than the typical soft scales of O(Aqcd), such as the 
normalization parameters of higher-twist DA's . The twist expansion is then "protected" 
by the combinations of l/m c and 1/r suppression factors, as one can explicitly prove 
by expanding LCSR in powers of l/m c (for more details on heavy-quark expansion of 
LCSR see [TIJ [23 EJJ ) . 

Another important step in the derivation of LCSR is the subtraction of the integral 
over higher states in This procedure introduces the lower limit uq = (m 2 —q 2 )/(sQ — 
q 2 ) in the convolution integrals entering LCSR (see the expressions for F^' t (q 2 , M 2 , Sq) 
in Appendix B). In LCSR's for D — > ir, K form factors the integration limit is no — 0.3 — 
0.4. Hence, the suppression of higher twists is not influenced, because after subtraction, 
the large part of the u-region in the integrals over DA's is retained. This is not the case 
in the limit of infinitely heavy quark and small q 2 , when, according to the scaling rule, 
Sq ~ (m 2 + 2m c u>) and uq ~ 1 — 2u>/m c — > 1. In this limit the power suppression of the 
higher-twist terms in LCSR depends also on the end-point behavior of DA's, modifying 
the initial twist hierarchy of OPE. For example, the formal l/m c suppression of the 
twist-3 terms mentioned above, is removed at m c — > oo. Still, this modification does not 
influence the numerical suppression of t > 4 contributions, provided the effective Borel 
scale r is kept large. 

The accuracy of the quark-hadron duality approximation is difficult to estimate in a 
model-independent way. In LCSR, one minimizes the sensitivity to this approximation, 
using not too large M 2 , in order to suppress the integral over the higher states expo- 
nentially. In addition, the effective threshold Sq is determined, calculating the L>-meson 
mass from the differentiated LCSR (|10p and adjusting the result to the measured value, 
as it was done in [12} [13] . This procedure is more reliable than fixing the effective thresh- 
old from a stability of the sum rule with respect to the Borel parameter variation. On 
the other hand, if the result of LCSR calculation actually reveals a weak dependence on 
M 2 (within the adopted interval), that is an important indication of the reliability of 
the method. 

3. Choice of the input 

In this section we specify and discuss the choice of input parameters entering LCSR's. As 
already mentioned, we use the MS scheme for the quark masses. In previous analyses, 
(e.g., in |15]I16|) the c-quark pole mass was used in the final sum rule, which is certainly 
less convenient for a correlation function with a virtual c quark. For the c-quark mass 
value we adopt the interval obtained from charmonium sum rules with 0(al) accuracy 

[IT], 

m c (m c ) = (1.29 ± 0.03) GeV, (12) 

2 The LCSR expressions in Appendix B are given in a compact form, using derivatives of DA's. To 
restore the inverse powers of M 2 one has to perform a partial integration. 
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where, conservatively, we double the error. This interval is in a good agreement with 
the recent lattice determination [25]. For the s-quark mass we take the interval 



m s {fx = 2 GeV) = (98 ± 16)MeV, (13) 

which covers the range of the QCD sum rule determinations with 0{o£.) accuracy |26j . 
Employing the well known ChPT relations |27| : 

„ 2m H „ 9 m\ — (m,, + m^) 2 /4 . , , 

R= — = 24.4 ±1.5, Q 2 = — , , ' = (22.7 ±0.8) 2 , (14) 

m u + m d m z d - mi 

the u, (i-quark masses and their sum can be calculated. Since we neglect m u and mj 
everywhere, except in the parameters fj, n and hk, we simply use the above relations and 
obtain (adding the errors in quadrature): 

<"< 2GeV > - 2,„,'("2'GeV) =(2 43±a42) GeV ' 

2 

H K (2 GeV) = mR ■ ^ = (2.42 ±0.39) GeV, (15) 



m s {2 GeV) 



i , i h r 2 -i 



with a remarkably small SU(3)fi violation. 

In our numerical analysis, the two-loop running for QCD coupling is used, with 
Oi s (mz) = 0.1176±0.002 [28] . The scale-dependence of the quark masses and parameters 
of DA's is taken into account in one-loop approximation. Furthermore, we use a uniform 
scale fi for all renormalizable parameters, with the same default value fj, = 1.4 GeV, as 

in previous analyses of LCSR's; note that parametrically, \i ~ <J m 2 D — m 2 . 

The hadronic inputs in LCSR's include the hadron masses m^o = 1.865 GeV, m n ± = 
139.6 MeV and m^± = 493.7 MeV [28J. The pion and kaon decay constants, f n = 
130.4 MeV and fx = 155.5 MeV [2H] normalize the twist-2 pion and kaon DA's, re- 
spectively. All other parameters of twist 2,3,4 DA's relevant for our calculation are 
collected in Table [S] in Appendix [XJ Let us briefly comment on our choice. Ex- 
pressing the twist-2 DA's in terms of Gegenbauer polynomials, we adopt the intervals 
af(l GeV) = 0.16 ± 0.01, aJ(lGeV) = 0.04 ± 0.01, determined in [13], by fitting the 
form factor calculated from LCSR to the measured shape of B — > it form factor. These 
intervals are in agreement with the other determinations of a^4 summarized in |19j . 
For the first Gegenbauer moment of the kaon DA we use = 0.10 ± 0.04, obtained 
in [18] from the two-point QCD sum rules with NNLO accuracy. Finally, the interval 
airf = 0.25 ± 0.15 is adopted [291 [19]. All other Gegenbauer moments are put to zero, 
the same approximation as in the previous LCSR analyses. In fact, as already noticed 
in [16], the sensitivity of the LCSR for the D — > tt form factor to Gegenbauer moments 
is less than in the B — > it case. In particular, we have checked that nonvanishing values 
of the next Gegenbauer moments and af, at the level of af and aj, respectively, 
have a small influence on the numerical results. For the twist-3,4 pion and kaon DA's, in 
addition to the normalization parameters already specified in (|15j) the set of parameters 
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determined and updated in [19] is used. Note that due to the smallness of f$ n ~ fox, 
the size of nonasymptotic corrections to the twist-3 DA's <j>^° ' ^ K is small. Hence, it 
is justified to take asymptotic DA's in the NLO twist-3 terms of LCSR calculated in 

pang. 

The remaining hadronic input in LCSR (jlOp is the decay constant fo- In previous 
applications of LCSR's to heavy-light form factors, e.g., in [J2J [JJH Q3], the two-point 
QCD sum rule for the heavy-meson decay constant was substituted in LCSR, leading 
to a partial cancellation between radiative gluon corrections. However, the two-point 
sum rule, with its own Borel-parameter range and effective duality threshold, introduces 
an additional uncertainty in the calculated form factors. Note that in the sum rules 
with MS heavy-quark mass the a s -corrections are not sizeable and, therefore, their 
partial cancellation is not that important. On the other hand, the decay constant fjj 
has already been measured in D —* Ivi [30] with a very good accuracy. For that reason, 
in our numerical analysis we prefer to use the experimental result, assuming the isospin 
symmetry: f D o = /d±, and taking f D + = 205.8 db 8.9 MeV from [30j . where we add the 
errors in quadrature. Importantly, this value is obtained assuming \V c d\ = \V US \, with 
\V US \ = 0.2255 db 0.0019 from [2B]. Extracting \V CC [\ below, we will take this into account. 
The two-point QCD sum rule prediction for fu used in previous analyses (e.g., in |15j ) 
agrees with the experimental interval, but has a larger uncertainty. 

Finally, we specify the "internal" parameters of LCSR (jlOf) : the interval of the Borel 
parameter M 2 and the effective threshold Sq . For the former, we choose the region 
M 2 = (4.5 ± 1.0) GeV 2 , close to the one used in [15] . The threshold parameter — 
(7.0 dz 0.5) GeV 2 is fixed by reproducing (within 2% accuracy) the .D-meson mass from 
the auxiliary sum rule, obtained from differentiating the LCSR (jlOf) over 1/M 2 and 
dividing the result by (|10p . With our choice of M 2 and Sq , the usual criteria are 
fulfilled for the LCSR ()10p : smallness of the subleading twist-4 contributions (< 5% of 
the twist- 2 term) and, simultaneously, suppression of higher state contributions (< 10% 
of the total correlation function). Since below we also calculate the form factors at 
negative q 2 , we checked that the adopted ranges of M 2 and Sq are equally applicable 
for -2 GeV 2 < q 2 < 0. For the second LCSR for (f+ n + f D J, the same intervals are 
taken for consistency. The differences between the Borel and threshold parameters for 
the sum rules for D —* ir and D — ► K form factors turn out to be negligible. Note also 
that the effective Borel scale r = M 2 /(2m c ) ~ 1.7 GeV is sufficiently large. 

4. Form factors at q 2 = 

Substituting in LCSR (fTUj) the input specified above, we calculate the form factors 
/o_(0) and /^-(O). The numerical evaluation was done in two different ways: firstly by 
a direct integration over imaginary parts of hard-scattering amplitudes, and secondly, 
applying the numerically equivalent method of analytical continuation explained and 
used in [151 122]. 

The results at the central values of input parameters are displayed in Table [TJ Their 
dependence on Borel parameter is shown in Fig. [TJ and exhibits a remarkable stability, 
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Figure 1: The Borel-parameter dependence of the form factors /^.(O) and /^^-(O) cal- 
culated from LCSR. 

even beyond the adopted interval. The scale-dependence displayed in Fig. [21 is also 
mild. Conservatively, we consider the variation of the calculated form factors with the 
scale change in the interval 1.0 < fi < 3.0 GeV as one of the uncertainties. In addition, 
we investigate the numerical hierarchy of various contributions to LCSR. The sample 
of results for fp^q 2 ) and fpxiq 2 ) is collected in Table [2j The dominance of the twist- 
3 LO contribution was anticipated, due to the factor fi^/nic > 1. At the same time, 
the subleading twist-4 contributions are numerically strongly suppressed. The NLO 
corrections to twist-2,3 terms are also small, a clear indication that the "soft-overlap" 
mechanism dominates in D — > tt, K form factors. 

Furthermore, we estimate separate uncertainties of our calculation by varying each 
input parameter within its allowed interval. All significant uncertainties of ^(0) 



Table 1: Form factors f^ n (0) and f^^iO) calculated from LCSR $10\) and the estimated 
uncertainties due to the variation of the input. 



Formf. 
centr. value 


M 2 




s 


(fD^exp 


m c 


m s 


Gegenbauer 
moments 


tw.5 
(est.) 


0.667 


+0.003 
-0.001 


+0.04 
-0.003 


±0.01 


±0.03 


+0.005 
-0.006 


+0.08 
-0.06 


±0.001 (aj j4 ) 


±0.017 


0.754 


+0.001 
-0.0004 


+0.04 

-0.006 


±0.01 


±0.03 


+0.005 
-0.007 


+0.09 
-0.06 


±0.003 (of) 
±0.03 (af ) 
±0.01 (a£ 4 ) 


±0.001 



0.8 
0.7 
0.6 
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Figure 2: The scale- dependence of the form factors fp n (0) and fp K (0) calculated from 
LCSR 



Table 2: Relative contributions to LCSR for the form factors f^iq 2 ) and f£)x(° 2 ) 



Contribution 


/5,(o) 


flM 1 = - 2 GeV^) 


/S*(o) 


fZ K tf = -2 GeV^) 


tw2 LO 


35.9% 


32.9% 


36.2% 


33.7% 


tw2 NLO 


6.3% 


8.4% 


6.0% 


7.9% 


tw3 LO 


66.0% 


59.5% 


67.9% 


58.9% 


tw3 NLO 


-9.5% 


-2.9% 


-10.1% 


-2.9% 


tw4 LO 


1.4% 


2.2% 


-0.07% 


2.4% 



are collected in Table [TJ (For brevity we do not show a similar table for /^/^(O) + 

fDn(K)(ty> presenting only the total uncertainty below.) The remaining small effects, 
e.g., due to the variation of a s , are not shown, but included in the total uncertainty. 
Note that, according to (fT5|) . the error related to both and hk is influenced by the 
uncertainties in the determination of m s and ChPT parameters. To roughly estimate 
the effects of the unknown twist-5 contributions to LCSR's, we assume that their ratio 
to the twist-3 contribution is equal to the calculated ratio of twist-4 and twist- 2 (LO) 
terms. Since in the sum rules for f^^ + f^^ there is no LO twist- 2 contribution, we 
conservatively assume that the magnitudes of the twist-5 and twist-4 terms are the same. 
In addition, to assess the effect of higher Gegenbauer moments a% and af on the D — > K 
form factor, we recalculated this form factor assuming af = ±of^ and = iaj. The 
small variations due to the abovementioned effects are treated as separate uncertainties 
and included in the error budget in Tabled) 
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Adding all uncertainties in quadrature, we obtain: 



= 0.671°;™ (16) 

f+ K (0) = 0.75±g:SS • (17) 

The second LCSR yields numerical results for the sums of the form factors: 

/^(0)+/^(0)= 0.461^, (18) 

&(0) + &(0) = 0.60l°;& (19) 

We also quote our prediction for the products of D decay constant and the form factors: 

/d/^(0) = 1371$ MeV, (20) 

/d/d*(0) = 155 i?5 M eV. (21) 



These quantities are independent of the experimental value of fr> , used to calculate (116]) 
and (|17p . and therefore have a slightly smaller uncertainty. 
Finally, the predicted ratio of the form factors is 

f+ (o) 

±I^-L = 0.88 ± 0.05 , (22) 

where the /o dependence drops out and some uncertainties largely cancel. 

In Table [3] we compare theoretical predictions for the form factors /^(O) and fp K {0), 
obtained in lattice QCD and from LCSR's. Our results are in a good agreement with 
the lattice determinations. The form factors (I16p and (|17j) and their ratio (I22p are also 
in accordance with the previous LCSR estimates [15, 16J. The D ^ K form factor is 
now more accurately determined than in [15], due to a better knowledge of the c- and 
s-quark masses and of the parameter of" . 

In the final part of this section we present our predictions for the slopes and ratios 
of the form factors at q 2 = 0, that is, at large recoil of the pion or kaon, (see 



Table 3: Comparison of theoretical predictions for the form factors f^ n (0) and fp K (0)- 



Method [Ref.] 


/&r(0) 




Lattice QCD [H] 

[32J 
[33J 


0.57 ±0.06 ±0.02 
0.64 ±0.03 ±0.06 
0.74 ± 0.06 ± 0.04 


0.66 ±0.04 ±0.01 
0.73 ±0.03 ±0.07 
0.78 ± 0.05 ± 0.04 


LCSR [15J 

m 

this work 


0.65 ±0.11 
0.63 ±0.11 

67+ - 10 

U - D ' -0.07 


0.75 ±0.12 
u - ' °-0.08 
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for definition). We start with the parameter 5 for D 
calculated by taking the ratio of two LCSR's: 



7r transitions, which is simply 



1 + 



0.69 ±0.09. 



(23) 



The slope of the scalar form factor at q 2 = normalized by f^ n (0) is another interesting 
characteristics of the form factor. We obtain: 



i i 



dq 2 



1.4 ±0.3. 



Combining these two parameters, we are able to predict the combination 

1 + l//3 D7r - <W = 1.02 ±0.18, 



(24) 



(25) 



consistent with the CLEO measurement 0:1 + I/Pd-k - $Dn = 0.93 ± 0.09 ± 0.01. The 
corresponding slope parameters for D — > K form factors predicted from LCSR are: 



S DK = 0.79 ±0.07, j3 DK = 1.6 ± 0.4 , 



(26) 
(27) 



so that 

1 + 1/Pdk-Sdk = 0.84 + 0.17 

also agrees with the experimental result [2]: 1 + 1/Pdk — $dk = 0.89 ± 0.04 ± 0.01. 

As discussed in the literature (see e.g, |34J), the values of 5 and (3 for heavy-light 
form factors reflect the proportion of their hard-scattering and soft-recoil components 
and, respectively, their deviation from the scaling behavior predicted in the combined 
heavy-quark and large-recoil limit. We postpone a more detailed discussion of these 
parameters to a future study. 



5. Determination of \V c d\ and \V CS \ 

The latest CLEO measurements of semileptonic charm decays [2], fitted to various form 
factor parameterizations (to be discussed in detail in the next section), yield the products 
fDn(0)\Vcd\ an d fDK(0)\V cs \. Both CKM matrix elements can now be determined using 
our predictions for the form factors /dtt(0) and /dx(0). 

Having an accurate experimental value for fo at our disposal, allows us to make 
the extraction of \V C( i\ less dependent on the theoretical uncertainty of LCSR, than in 
previous analyses, where a sum rule prediction for fo was used. We employ the CLEO 
result [30J for the D-meson decay constant multiplied by |V^|> (i-e, without the additional 
assumption \V cd \ = \V US \): 

f D \V cd \ = 46.4 ± 2.0 MeV. (28) 

On the other hand, the CLEO result [2] for the same CKM matrix element multiplied 
by the form factor is : 

/D7r(0)|V cd | = 0.150 ± 0.004 ± 0.001 , (29) 
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obtained from the fit of the g 2 -bins in D — ► -Rev in a form of the series parameteriza- 
tion. The product of the above two experimental numbers is then divided by the LCSR 
prediction (|20p . yielding \V c d\ 2 , from which we obtain: 

\V cd \ = 0.225 ± 0.005 ± 0.003 to™ , (30) 



where the first and second errors originate from the experimental errors in (|28|) and 
(|29|) . respectively, and the third error is due to the uncertainty of LCSR. Note that this 
procedure involves the square of |V^| on the experimental side. Hence, the theoretical 
uncertainty of LCSR given in (I20p approximately halves in (|30p . Our result is in a good 
agreement with the value \V c a\ = 0.234±0.007±0.002±0.025 determined in [2] by using 
the lattice QCD value of f^ )n (0) from [32J. This agreement is not surprising because the 
form factor obtained from LCSR is close to the lattice result (see Table [3|). 
Furthermore, we determine the ratio of \V c d\ to | V^ s | , dividing (j29j) by 

fDK(P)\Ves\ = 0.719 ± 0.006 ± 0.005, (31) 

obtained from D — > Kev e data fit [2j. Using the ratio (|22[) we obtain: 

j^j = 0.236 ± 0.006 ± 0.003 ± 0.013 , (32) 

I 'cs I 

where the first and second uncertainties are due to the combined (in quadratures) errors 
in (|29|) and (I3ip . respectively, and the third uncertainty stems from the LCSR calculation. 
Within errors, this ratio is in agreement with |V^|/|T4 S | obtained from the values quoted 
in [2], where |V^ S | = 0.985 ± 0.009 ± 0.006 ± 0.103 was determined using the form factor 
/ojf(0) from lattice QCD [32 j . Our determinations (|30p and (|32p are consistent with 
\Vcd\ = \V US \ and \V CS \ = \V ud \. 



6. Form factors and their shapes at q 2 ^ 

The form factors /^^(q 2 ) and f^ K (q 2 ) are analytic functions of the complex variable q 2 . 
The singularities located on the real positive axis include the poles of the ground-state 
vector mesons D* and D*, respectively, and their radially excited states. In addition, 
there are branch points, generated by the hadronic continuum states, starting from Dtt 
and DK thresholds at q 2 = (mp +m n ) 2 and q 2 = (mjj +m^) 2 , respectively. (Note that 
the D s ir intermediate state in the D — > K form factor is forbidden in the isospin limit). 
To obtain convenient parameterizations for the form factors, one employs analyticity in 
two different ways. 

The first approach uses the dispersion relation, e.g., for the D — > ir form factor: 

oo 

f + D M 2 ) = : % 2 + / ds?p^, (33) 
1 — q jra D * J s — q A 

(rnD+niTr) 2 

where the ground-state D*-pole at m 2 -,* = (2.01 GeV) 2 [2H] is isolated, and the hadronic 
spectral density p^ n (s) includes all other intermediate hadronic states with the D* 



15 



quantum numbers. In fact, fn 2 Dt is slightly larger than the continuum threshold (so that 
the D* — > Dtt decay is observable), whereas in the dispersion relation for ffixiq 2 ) ^ ne 
pole at rnj-ft = (2.112 GeV) 2 [28J lies below the threshold. Importantly, the dispersion 
relation (|33j) has no subtractions, due to the expected QCD asymptotics of the form 
factor lim| (J 2|_^ 00 f^) n (q 2 ) ~ 1/|<? 2 |- The residue of the pole in (|33|) is normalized as: 

f D *g D * Dn 

CD * = IST ' (34) 

where fjj* and go* Dw are the D* decay constant and D*Dir strong coupling, respectively, 
defined in the standard way (see e.g., [9]). 

The rigorous dispersion relation ([33]) is valid at any q 2 . Hence, matching a calculated 
form factor, e.g., the one obtained from LCSR, to the dispersion relation in the region 
where this calculation is valid, one can, in principle, predict the form factor outside this 
region. This is, however, only possible, if the complicated integral over the hadronic 
spectral density in (|33p is parameterized in a simple and reliable way. One possibility is 
to replace this integral by an effective pole: 

JdM 2 ) = °d* f 1 TT^ 1 TFT 1 2~t) ' ( 35 ) 

where aon and jd-k parameterize the residue and the position of this pole, so that the 
normalization of the form factor at q 2 = is: 

f+ w (0) =c D *(l-a Dn ). (36) 

Using the relation 7^ = 1/aon inspired by the combined heavy-quark and large-recoil 
limit, the two-pole ansatz is reduced [20J to the specific BK-parameterization: 



*DkU ) (1 _ q 2/ m 2 Df){l _ a Dv q 2 /m 2 D ,) " {6 ' } 

The second approach based on the analyticity of the form factors employs the confor- 
mal mapping of the g 2 -plane (see e.g., |35| for the early uses): 



^to) = ^pSH== ^ 
\A+ - q + Vt+ - *o 

where t± = (mp ± m^^^) 2 and to < t + is an auxiliary parameter. Applying this 
transformation, one maps the g 2 -plane (with a cut along the positive axis) onto the inner 
part of the unit circle in the z plane, so that fD-K(K)(q 2 )q 2 ^z is f ree from singularities 
at \z\ < 1. The conformal mapping (|38|) was employed while deriving the unitarity 
bounds for the heavy-light form factors in |36] I37j. Independent of these bounds, with 
an optimal choice of to, the semileptonic region < q 2 < (mp — rn n (K)) 2 i s mapped 
onto the interval of small \z\. Hence, a simple expansion in powers of z around z = 0, 
retaining only a few first terms, provides a reasonably accurate parameterization of the 
form factor. 
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In this paper we will use the recently suggested version [21] of the series parame- 
terization for B — > 7r form factor. Adapting it to the case of D — > tt transition, we 
have: 

1 N ~ 

fLil 2 ) = i a y m2 i^Mf- (39) 

1 1 l m D* k=0 

As explained in [21], this parameterization ensures general analytic properties of the 
form factor: the D*-pole, branch point at q 2 = t + and ~ 1/q 2 asymptotics at large q 2 . 
Furthermore, to obey the expected near-threshold behaviour, the relation 



b N 



( ^\N N -' L 



(40) 



k=0 



has to be introduced, reducing the number of independent parameters by one. In ad- 
dition, we find it more convenient to keep the form factor at zero momentum transfer 
/d„.(0) as one of the independent parameters, correspondingly rescaling the coefficients 
in the series expansion: = f^j n (0)bk, so that 



N-l 



bo = 1- ^2 h 



fc=i 



k 



z(0,t Q r -z(0,t ) 



N 



(41) 



This leads to the final form of the series parameterization used in our analysis: 

N-l / 



1 — q 2 /m 2 Dt 



l+Y,bk[z(q 2 ,t ) k -z(0,t ) k 



k=l 



k-N 



N 



z(q 2 ,t ) N -z(0,t ) 



N 



(42) 



In [21] the advantages of this choice with respect to the previous versions [36] are dis- 
cussed (see also [3H]). The analogous ansatz for f^ K contains m 2 Dt in the pole prefactor. 
Note that the formal prescription for the conformal mapping is to multiply the r.h.s. 
of (|42|) by the pole factor if the ground-state pole lies below the threshold t + , which 
is, strictly speaking, only valid for D — > K form factor. We prefer to retain this factor 
in (|42p also for D — * n case, having in mind that D* is located very close to the Dtt 
threshold. In what follows, we simply rely on the smallness of the variable z in (I42p . not 
taking into account the unitarity bounds [36] for the coefficients b^, because, at least for 
small N, these bounds are not restrictive [21] (see also [37]). 

In order to match the LCSR prediction for D — > tt and D — > K form factors to 
one of the parameterizations discussed above, we have to calculate these form factors 
beyond q 2 = 0. However, the small part < q 2 <C m 2 of the semileptonic region 
where our calculation is valid, is too narrow to serve as a "lever arm" for fitting various 
parameterizations. 

In this work, in order to enlarge the interval of q 2 , we calculate the form factors at 
Imin < q 2 < 0, that is, at negative momentum transfers not accessible in semileptonic 
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decays □. Fitting the LCSR predictions in this region to a certain parameterization, we 
then use the analyticity of the form factors, continuing the fitted parameterization to 
positive q 2 and accessing the whole semileptonic region < q 2 < t_. 

Note that LCSR's are fully applicable at q 2 < 0, since the virtuality of the c quark 
in the correlation function is even larger at q 2 < 0, than at q 2 = 0. For our purpose 
it is sufficient to take q 2 not too large, and the actual numerical calculation is done 
up to q 2 min = —2 GeV 2 . In fact, there are several reasons to keep moderate values of 
\q 2 \. First, we still can use the same ranges of the sum rule parameters M 2 , s®, fi, as 
specified in Sect. 3, whereas at very large \q 2 \ some of the choices have to be modified. In 
addition, at large virtualities, \q 2 \ 3> « 2 ,m 2 , large logarithms in NLO terms can destroy 
the balance of perturbative expansion. Finally, at very large negative q 2 , the lower limit 
of integration no in LCSR moves too close to 1, and this may potentially influence the 
twist expansion. 

Turning to the fit procedure, we fix the q 2 = values of the form factors (|16p and 
(|17p obtained from LCSR's and, in addition, calculate the shapes of the form factors 
fDir(K)(l 2S } '/ ' ^d-k(k)^ a ^ — ^ GeV 2 < q 2 < 0. The uncertainty is determined in the same 
way as described in Sect. 4 for the form factors at q 2 = 0. The calculated D — ► it 
and D — > K form factors and their shapes at q 2 < are displayed in Figs. 3 and 
4, respectively. Note that the shapes have smaller uncertainties than the form factors 
themselves, because in the ratios some uncertainties cancel (e.g., the one due to fjj). 

The calculated shapes, taking into account their uncertainties, are then fitted to the 
BK and series parameterizations, presented in (|37p and (|42p . respectively. In the series 

parameterization (j4~2j) we choose to = t + — ^/t + — t_yjt + — q^~ n -, so that the whole inter- 
val q^in = —2 GeV 2 < q 2 <t- relevant for the calculation and subsequent continuation 
of the form factor is mapped onto the narrowest possible interval \z\ < 0.22 (\z\ < 0.09) 
for D — > 7r (D — > K). We found that both parameterizations describe the shapes cal- 
culated at q 2 < reasonably well. In addition, we also fitted the form factors to the 
simplest one-pole parameterization, but the fits yield an unnaturally small, lower than 
rriD* , pole mass. Moreover, after continuing to positive q 2 , the one-pole form factor 
noticeably deviates from the experimentally measured shape. 

The parameters of BK parameterization old-k and oldk obtained from our fit are 
shown in Table HI in comparison with the experimental [2] and lattice QCD [32J BK-fits. 
Previous LCSR estimates of these parameters |15j are smaller but have also larger errors; 



Table 4: The shape parameters of BK parametrization 



Method 


Ref 






LCSR at q 2 < 


this work 


n91 +U.ll 


OiDK — U.17_ 13 


Experiment 


m 


a Dw = 0.21 ±0.07 ±0.02 


a DK = 0.30 ±0.03 ±0.01 


Lattice QCD 


L32J 


a Dn = 0.44 ±0.04 ±0.07 


a DK = 0.50 ±0.04 ±0.07 



3 In fact, this region corresponds to a hypothetical (but still physical) process of Ity(K) — > viD scattering. 
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Figure 3: Form factors /^(q 2 ) (upper panel) and f^ K (q 2 ) (lower panel). The LCSR 
results with uncertainties (points with error bars at q 2 < 0) are fitted to 
series-parameterization (solid line and dashed lines indicating uncertainties) 
and compared to the CLEO measurements J7}/ (points with error bars at q 2 > 0), 
where \V c d\ and \ V CS \ are taken from [28]. 



in fact, they have been determined from a different procedure, where, in addition to the 
LCSR form factor at q 2 = 0, the calculated D*Dir coupling was used, adding its own 
uncertainty. 

For the series parameterization (I42|) . a fit is possible already at N = 2, i.e., with only 
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one free parameter b\ for the shape (62 is fixed from the condition (|40p ). We obtain: 



6f* = -0.8ljj3 , bf K = -0.9±°; 8 7 . (43) 

Fits at A" = 3, 4 were also performed, yielding numerically very close results. For N > 5 
the unitarity bounds [21] start to constrain the coefficients bk- 

In what follows, we choose the N = 2 series parameterization to be our preferred 
analytic expression for the shape, having in mind that it is less model-dependent than 
the effective pole ansatz for the dispersion integral. Continuing this parameterization to 
the semileptonic region q 2 < t- = {tub — i^n(K)) 2 = 2.98 GeV 2 (1.88 GeV 2 ), we compare 
in Fig. [3] the form factors with the experimentally measured ones, presented in pQ in 
q^-bins. For the normalization of the data we take the averages: \V c d\ = 0.230 ± 0.011 
and \V CS \ = 1.04 ± 0.06 from [28J. The form factor shapes, which are independent of 
normalization at q 2 = and CKM parameters, are displayed in Fig. [H We compare our 
predictions for the series parameterization with the shapes obtained in [2] and observe 
a good agreement 0. 

Furthermore, we calculate the total semileptonic widths divided by the square of CKM 
parameters from 



(mo-m^) 2 

W^\ 2 ~ 24^3 



2m p 



\HM 2 )W (44) 







(at mi = 0) and the analogous formula for Y(D —>■ K~£ + V£)/\V CS \ 2 , using the predicted 
shape of the form factors (with the series parameterization) and their normalization at 
q 2 = 0. Again, in the case of D — » tt a better accuracy is achieved by normalizing with 
the product /^/^(O) calculated from LCSR. To this end, multiplying both sides of (|4"4"|) 
by f D (in the isospin limit) we replace this factor on l.h.s. by the leptonic width, using 

r(D+ -> £+u e ) = ^k\ Vcd \ 2 f 2 D m 2 m D ( 1 - V , (45) 



07T \ m D 



and obtain the following prediction: 



(4.7lJ;|) • 10" 28 Gey 2 . (46) 



Employing the experimental numbers for the branching fractions from the latest CLEO 
measurements: BR(D° -+ Tr~e + u e ) = 0.288 ± 0.008 ± 0.003% [2], £i?(L>+ -> = 
(3.82 ± 0.32 ± 0.09) x 10" 4 [30], and using t d ± = (1.040 ± 0.007) ps, t d o = (0.4101 ± 
0.0015) ps [28], we obtain 



= 0.221 ± 0.002 ± 0.005 ^Zgloxi . ( 47 ) 



4 Note that in [5] a different version of series parameterization is used to fit the shape, hence we do not 
directly compare the fitted parameters. 
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Figure 4: The shapes of the D — > 7r (upper panel) and D — > if (lower panel) form factors 
obtained from LCSR's at negative q 2 (points with error bars), fitted to the 
series parameterization (solid line and dashed lines indicating uncertainties) 
and compared to the shapes measured by CLEO fMj (shaded regions). 



with the errors originating from the semileptonic and leptonic branching fractions and 
theoretical uncertainty, respectively. This determination is somewhat independent of 
([30]) because it involves also the measured shape of f + (q 2 )\V c d\- For the ratio of total 
semileptonic widths we obtain: 

IVU 2 T(D° -> ir-£ + v e ) , . 

1 1.65 ±0.2. (48) 



\V cd \ 2 T(D^K-l+u,) 
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Substituting the CLEO results for the branching fractions of these channels: the one 
quoted above and BR(D° -» K~e + vi>) = 3.50 ± 0.03 ± 0.04% [2], yields 

= 0.223 ± 0.003 ± 0.002 ± 0.015 . (49) 

where the errors are from semileptonic D — > ir, D — > K branching fractions and the 
LCSR result ([Ml), respectively. To shorten this discussion, we do not present here a 
comparison with the data of other experiments on charm semileptonic decays [H EJ [6] , 
having in mind their general agreement with the CLEO data. 

Concluding this section, we turn to the scalar form factor /xw^nC? 2 ) which is obtained, 

substituting the LCSR results for f£ <K) {q 2 ) and [f^ {K) (q 2 ) + f^(K)(l 2 )} in ©• We 
fit the scalar form factors calculated at negative q 2 to the series parameterization of the 
type (|42p (without the D^-pole factor which is irrelevant in this case). The results for 
N = 2 are 

b{°^ = -2.6l°;f, b{°> DK = -3.3l°j . (50) 

Going to N > 3 demands dedicated unitarity bounds in the scalar heavy-light channel 
which, to our knowledge have not been derived and are beyond our scope. 
The predicted scalar form factors are plotted in Fig. 

Our results are in agreement with and fp K presented in [32] in a form of BK 
parameterization 

fDTr(K)(l 2 ) = ~-\ TTTa 2 \ ' 

1 > l-q 2 /{P D7r ( K )m z D * ) 

with Dn = 1.41 ± 0.06 ± 0.07, (3 DK = 1.31 ± 0.07 ± 0.13. Our fit to (J5U) yields: 
Pd-tt = 1-25 ± 0.2, Pdk = l-3lo3, w ith larger uncertainties. 

Finally, we extrapolate the scalar form factors to the unphysical point q 2 = m? D , 
located slightly above i_, and obtain 

fUml) = 1-40±°;?1 , f DK {m 2 D ) = 1.29^ . (52) 

These results can be compared to the approximate relation 

, lim , fL(K)(<l 2 ) = fD/L(K) = 1-58 ± 0.07 (1.32 ± 0.06) , (53) 

q 2 — >m A D 

derived from the current algebra combined with the soft pion (kaon) limit |39j . where 
we used the measured values [281 ED] of the decay constants. 



7. Summary 

In this paper we returned to the calculation of D — > ir and D — > K form factors from 
LCSR's. Several improvements have been implemented, including the use of MS c-quark 
mass, and updated parameters of pion and kaon DA's. 



22 




q 2 [GeV 2 ] 




-2.0 -1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 

q 2 [GeV 2 



Figure 5: Scalar form factors /j^on obtained from fitting the LCSR result at negative 
q 2 to the series parameterization. Notations the same as in Fig. d 



The main advancement is in the phenomenological direction. Employing the accurate 
measurement of the .D-meson decay constant, we effectively decreased the theoretical 
uncertainty of \V c d\ determination from the D — > ixlv\ decay distribution, using the 
LCSR prediction for the product of the D — > tt form factor and fo- The uncertainty 
in the determination of |V^ S | is also reduced, due to a better knowledge of the s-quark 
mass and various SU(3)fi violating effects in the kaon DA's. Our results for \V c d\ and 
| V cs | are in agreement with lattice QCD determinations. 

A new element presented in this paper is the prediction of the D — > tt, K form factors 
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in the whole semileptonic region, combining LCSR calculation with the analyticity of the 
form factors. The latter property is cast in the form of conformal mapping and series pa- 
rameterization, in the version recently suggested in |21j. The form factor shapes obtained 
from this combined procedure are in a good agreement with the latest experimental mea- 
surements of the semileptonic charm decay distributions by CLEO collaboration. Our 
analysis based on the conformal mapping can be further refined, by using more terms in 
the power series and implementing the constraints from the dedicated unitarity bounds. 
Applications to other hadronic form factors calculated from LCSR are also possible. 

Another interesting task which will be studied elsewhere, is the comparison of B —* n 
and D — > tt form factors, calculated at two different finite quark masses, with the com- 
bined heavy-quark mass and large recoil limit. Such a comparison will allow one to 
quantify the deviations from the symmetry relations in this limit, as well as the propor- 
tion of hard-scattering and soft-overlap mechanisms in the heavy- light form factors. 

LCSR's provide analytical, but essentially approximate expressions for the hadronic 
form factors. The accuracy of light-cone OPE is limited, due to finite amount of terms in 
the twist expansion and uncertainties in the parameters of pion and kaon DA's. A further 
improvement of OPE is possible, e.g., if the subleading twist-5 terms are calculated. For 
that one needs a separate study of pion and kaon twist-5 DA's. The gluon radiative 
corrections to the subleading twist-4 and three-particle terms represent a technically 
difficult task, but we expect these corrections to be very small. A further limitation of 
the accuracy is caused by the quark-hadron duality approximation used to model the 
hadronic spectral density in LCSR. The resulting uncertainty is difficult to estimate, 
still it is effectively minimized in the sum rules. 

The comparison of our predictions for charm semileptonic decays with experiment and 
lattice QCD ensures optimism and provides an additional test for the important applica- 
tions of the LCSR method, such as the \V U b\ determination from exclusive semileptonic 
B — > 7r decays. 
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Appendix 



A. Summary of the A -meson distribution amplitudes 

Here we present the definitions and the expressions of the two- and three-particle K- 
meson DA's of twist 2,3 and 4 used in LCSR. The corresponding formulae for pion DA's 
are obtained by replacing everywhere K — > it, s — > d, mx — ► rn n , m s — > — > and 

—> Hit- 



A.l. Definitions 

The K-meson two-particle DA's are defined from the following bilocal matrix element: 
(K-(p)\sUxi)ul(x 2 )\0) x ^ = l -^f K j\u e^+^([pj 5 ]^ K (u) 

+^[?*75]^0i - x 2 ) 2 (p4K(u) - - ^2)75]^ J ip4x(v)dv J , (54) 

' 

where the variable u (u = 1 — u) is the fraction of the meson momentum carried by the 
s-quark (light antiquark). This decomposition contains the twist- 2 DA (Pk(u>), twist-3 
DA's 4> p ^ K {u) and <j>%^{u), and twist-4 DA's 4>4k(u) and ^ak{u)- The definition of each 
separate DA is easily obtained, multiplying both parts of the above equation by the 
corresponding combinations of 7 matrices and taking Dirac and color traces. 
In the same way, the three-particle DA's are defined via the matrix element: 

(K-(p)\t(xi)g s G%(x 3 )ul(x 2 )\0) x2 ^ = ^ J V a i e i ^+ aM ^ 



~ e,vS P (lx)^{(p p 9 SX ~ P 5 9 pX )*Mm) + Pa( ^. X&) (*4K:(ai) + *«-(«<)) } 

(55) 



where $3E"(ai) is of twist-3 and the other four DA's are of twist-4. 

In addition, in [19] one more three-particle twist-4 DA (aj) is taken into consid- 
eration, originating from the operator which contains a covariant derivative of the gluon 
field: 
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(K(p) \S(x 1 ) 1 ^gD a G a p{x z )q{x 2 )\ 0) = -if K p^ J V Ui e ip(^x+^2+a 3 x 3 ) Em ( a .y 

(56) 

The three-particle DA's depend on momentum fraction variables 
oti = {ai,«2!«3} and the integration measure is Vcti = da.1da2da.35 {I — a\ — a 2 — 03). 

For the total antisymmetric tensor we use the convention e 0123 = — 1 5 which corre- 
sponds to Tr{-i b ^Yl a l P } = Aie^ 13 . 

In the following, we list the expressions for all DA's entering (|54|) . (|55|) and (I56p . based 
on NLO in conformal expansion and operator identities and updated in [19|. 

A. 2. Twist-2 distribution amplitudes 

The twist-2 (px(u) distribution amplitude is expanded in a series of Gegenbauer poly- 
nomials: 

<p K (u, ») = 6u(l-u) (l+ Yl «n(^)Cn /2 (2n-l) J , (57) 

V n=l,2,- J 

where only the first two coefficients (Gegenbauer moments) af (/u) and af" (p) are retained 
and LO scale-dependence is taken into account. The formulae for the scale-dependence 
of these and other relevant DA parameters can be found, e.g., in [19J. 

A. 3. Twist- 3 distribution amplitudes 

In the same approximation, the twist-3 DA's are described by px and three additional 
parameters f3K,^3K, <^3K- Their definitions in terms of hadronic matrix elements of 
local operators are given in |19j . We also use the short-hand notation 

Hk 

V3K = 7 • 

JKPK 

In our calculation we neglect the u, d quark masses, hence the expressions presented 
here are somewhat simpler than the original ones in [19) . In particular, in the adopted 
approximation the parameters p+,p^; introduced in [TS] are equal: 

p« = p K _= P K = ^. (58) 

PK 

The twist-3 kaon DA's used in our calculation are: 

f ZK (u) = l + 3p K (l-3af + 6af)(l + lnu) 

K 

- (3 - 27af + 54af ) C{' 2 (2u - 1) 

+ 3^Wrj 3K -p K {af -5af))c 2 1/2 (2u-l) 

+ (lOvsxhK ~ \p K 4^j Cl /2 (2u - 1) - 3 mK uj 3K Cl /2 (2u - 1) , (59) 
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4% K (u) = 6uu< 1 + !-— (3 - 15af + 30af ) 



15 



+ p *(3af-^af )cl'\2v- I) 



+ 5 f W(10 - uj 3K ) + 3p K aA Cl ,2 {2u - 1) + ri 3K X 3K cl /2 (2u - 1) 
+ 3p K (1 - 3af + 6af ) In u I , (60) 



®3K(ai) = 360aia 2 a z 3 <! 1 + A 3 ^(ai - a 2 ) + uj 3K - (7a 3 - 3) 



(61) 



A. 4. Twist-4 distribution amplitudes 

As explained in detail in |19j . the twist-4 DA's are described by 13 parameters of 
the conformal expansion. They are expressed via three nonperturbative parameters 
S'jf , ui\k , k^k and in addition fixed by the renormalon model of twist-4 DA's. Here we 
give the expressions for the twist-4 kaon DA's, where the above mentioned relations are 
already substituted and the same approximation as for twist-3 DA's is adopted. 

We rederived the expressions (4.27) and (4.28) in |19j for the twist-4 two-particle DA's 
defined in (4.26) and using the operator relations given there in (Al),(A2). We found that 
both (4.27) and (4.28) should be corrected by replacing tpi-K(u) — * ^a-k{u) + mj^^K 
(in the notations of [IDJ). In fact, our version of i(^4-k(u) agrees with the function B(u) 
introduced in [3D]. Moreover, we restore the correct normalization ip^K(u)du = 0. 

We use the following expressions for two-particle twist-4 DA's: 



Mu)=1 pT4 {u)+xp WW {uh 



where 



^ 2 
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K 



Cl /2 (2u 



l) + f«fC 3 1/2 (2n-l) 



(62) 



(63) 



and (the corrected version) 



<K W (U) 



+ 



+ 



m 2 K < 



6^(1 -3af +6a£) 



K 



Cl ,2 {2u-l) 



18 

— af + Zp K {l - 9af + 18af ) + 12k ak 



C\ ,2 {2u - 1) 



2-6p K (af -5af)+60 mK 



K \ 



Cl /2 (2u-l) 



_l r 18 „K n„K „K 16 



9p a 2 + — k ak + 20r] 3K X 3K 



C 3 1/2 (2n-l) 



9 



C\ ,2 {2u- 1) ) +6m^(l-3af + 6a£)lnu, (64) 



K 
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Table 5: Parameters of the pion and kaon DA 's (normalized at 1 GeV): of, a\ are fitted 



in f73j/ ; a± from J7&]/ . all others from fig) /, calculated from {7 Sty 



a l 


o 


a? 


10 ± 04 


"2 


16 ± 01 




25 ± 1 5 


a\ 


0.04 ± 0.01 












" >4 





J3 


(0.0045 ± 0.0015) GeV^ 


/s* 


(0.0045 ±0.0015) GcV^ 


<4 


-1.5 ±0.7 




-1.2±0.7 


AJ 







1.6 ±0.4 




0.2 ±0.1 




0.2 ±0.1 




(0.18 ± 0.06) GeV 1 ' 


& 


(0.2 ± 0.06) GeV^ 









-0.12 ±0.01 



where 



^K(u) = <pl*( U ) + ^ K w (u 



WW, 



(65) 



6 2 



K 



/ 200 



+ 196(2u - l)a{ 



K ^ u 2 u 2 



u <— > u 



\u «-> u 



+ 21w 4A - ( im(2 + 13uu) + [2u 3 (6u 2 - 15u + 10) Inu] + [ 

- 14af ^uu(2u - 1)(2 - 3uu) - [2u 3 (u - 2) lnu] + [ 
= m 2 K l ^K iK (uu(2u - 1)(1 - 2uu) + [5(u - 2)u 3 lnu] - [ 
+ 4?73^uu^60n + 10A 3 ^ [(2u - 1)(1 - uu) - (1 - 5uu)] 

- uj 3K [3 - 21uu + 28uV + 3(2u - 1)(1 - 7uu)]\ 



(66) 
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K 



itit (4 - 9im ± lKkrtr) + [u d (10 - 15tt + 6it 2 ) lnu] + [u ^ u) 



K 



+ Auu(l + 3uu) ( l + -(2u-l)af 



The twist-4 three-particle DAs have the following expressions: 

f ^ 21 ' 7 .AT, 



(67) 



120aia 2 a3S -^-("l — "2)^4^ + — a{ (1 - 3a 3 ) 



20 



9 / , ]f 1 

— (ai - a 2 )a 2 + -k 4 k 



(68) 
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120aia 2 a3^|- + - af (ax - a 2 ) + —^4^(1 - 3a 3 )|, (69) 
1 



30a| I 81 



21 



-(«! - a 2 ) + —of [ - «3(1 - 03) + 3(ai - a 2 ) 2 ] 



+ — w 4j ftr(ai - a 2 )(l - 2a 3 ) 



+ m#(l - a 3 ) 



-K4K 



(70) 



*4Ar(«i 



21 



1 7 

-(1 - a 3 ) - —af(ai - a 2 )(4a 3 - 3) 



+ — w 4 #(l - a 3 )(l - 2a 3 ) 



+ m 



K 



^af(a? - 4aia 2 + a 2 ) - ^(ai - a 2 )K 4K 



(71) 



(72) 



Using the equations of motion, the parameter kak can be expressed via of and the 
quark mass: 

^ = -8"40 ai + W (73) 
and the parameter entering 'E^k is taken from the renormalon model |19j: 



-0 — gOE" ! 



(74) 



The numerical values for all parameters entering the pion- and kaon- DA's are collected 
in Table E 

B. Contributions to LCSR 

Here we present the separate contributions to LCSR's for the form factors f^ K and 
(fr>K + fix) including LO twist-2,-3 and -4 terms. The corresponding contributions 
for D — > 7r form factors are obtained by replacing: the kaon DA's by the pion DA's, 
rriK — ► Tn n ~ and m s — > rrid ~ 0. 



F«>\ q \M\s°) 



1 

J u 



2 2 - 2 - 
uM 2 



(75) 
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F K >\q 2 ,M 2 ,s°) 



m 2 J K 



du e 



2 2-, 2 
m c —q u-\-nij^uu 

uM 2 



UK 



UQ 




m 2 + q 2 



u m 



K 



d 



2um 2 K m 2 



2{m\ 



+ u 2 m? K ) du (m 2 — q 2 + u 2 m 2 K ) 2 



m: 



urn 



K 



+ 



w \ mi 



3m 2 K 



m 2 — q 2 + u 2 m 2 K 



q 2 + u 2 m z K 

d_ 

du 



d 

du 



2um 2 K 



m 2 — q 2 + u 2 m 2 K 




Izk{u) 



(76) 



_L 

F K > 4 (q 2 ,M 2 , S °)=m 2 J K j due 

UQ 



u+m jrUU 



m 2 — q 2 + u 2 m 2 K 



2u 2 m 2 K 



uiPak(u) + ( 1 9 "V7 K 9 2 ) / dvif) 4K {v) 

1 — q l + u z m z K J J 



mi 



m 2 u 



(*— 



&um 2 K 



d 



+ 



\2um K 



4(m 2 — q 2 + u 2 m 2 K ) \du 2 m 2 — q 2 + u 2 m 2 K du (m 2 — q 2 + u 2 m 2 K ) 2 
d_ 
du 



(«) 



2um 2 K 



m 2 — q 2 + u 2 m 2 K 



h K (u) 



2um\ 



m 2 — q 2 + u 2 m 2 K \ du 



2™2 



u— + 1 



du 
^u 2 m 2 K 



mi 



q A + i^m 2 ^ 
Gwb 2 ^ 



2um 2 K {m 2 — q 2 



u m 



K) 



(m 2 — q 2 + u 2 m 2 K 



2 \2 



/4a:(m) 



cfot m 2 — g 2 + u 2 m\ 



K 



— 12 

™>lf K e ^ 



4(m 2 — g 2 + m 2 K ) 2 



du 



u— »1 



(77) 



where u = 1— u, iio = ^ 2 — + m lf + ^/ ( s o* ~ Q 2 ~ m \) 2 + ^m 2 K {m 2 — q 2 )j /{2m 2 K ), 
and the short-hand notations introduced for the integrals over three-particle DA's are: 

u 1 



= j dai j 



dOL\ J dv &3K(o!i 
(u-ai)/(l-ai) 



«2 = 1 — «1 — «3 5 

03 = (ti — a\)/v 



(78) 
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Hk{u) 



Iak(u) 



lf K (u) 



U 1 

(u-ai)/(l-ai) 
u 1 



'<• / < l <\ l / ^ (2u - 1) $3K(oti 



oli = 1 — ol\ — as, 

«3 = (u — U\)/v 



doti 

(u-ai)/(l-ai) 



(u-ai)/(l-ai) 



2*4i^(ai) - $4^(a*) + 



«2 = 1 — ai — «3, 

CK3 = (li — &l)/v 



«2 = 1 — «1 — «3 5 

«3 = (-u — a±)/v 



u 

J da.\ 

(u-ai)/(l-ai) 



u(l - v)E 4K (ai 



ct2 = 1 — ol\ — as, 
as = (u — a\)/v 



The LCSR for (fp K + f^x) i n LO has the following contributions: 

F^ 2 (q 2 ,M 2 ,m 2 c ,s^)=0, 



/2 2 — i 2 i 
m r .—q u+mj^uu I < , 
due' um'^ M 



+ 



\fKfnc) rn 2 — Q 2 + u 2 m? K \du m 2 — q 2 + u 2 m\ 



2um 2 K 



6u du 
hK{u) \ , 



where 



U 1 

h K {u) = Jd ai J ^ [(3 - 2v)] * 3A -(ai 



(tt-ai)/(l-ai) 



«2 = 1 — ol\ — as, 
as = (u — a\)/v 
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and 

i 

,2 



F^(q',M z ,s^)=m z c f K / due SS» 



2 2- , 2 - 
m c —q u-\-rrij^uu 



2um 2 K f u 2um 2 K 



dvijj4K(v) + 



+ u 2 m 2 K Jq m 2 — q 2 + u 2 m 2 K \du 2 



m c — 9 2 + u 2 m 2 K du (m 2 — q 2 + u 2 m 2 K ) 2 



(86) 
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